Electronic structure of nickelates: from 2D heterostructures to 3D bulk materials 
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Reduced dimensionality and strong electronic correlations, which are among the most important 
ingredients for cuprate-like high-Tc superconductivity, characterize also the physics of nickelate 
based heterostructures. Starting from the local density approximation we arrive at a simple two- 
band model for quasi two-dimensional (2D) LaNiOs/LaAlOa heterostructures and extend it by 
introducing an appropriate hopping in the z direction to describe the dimensional crossover to three 
dimensions (3D). Using dynamical mean field theory, we study the effects of electronic correlations 
with increasing interaction strength along the crossover from 2D to 3D. Qualitatively, the effects 
of electronic correlations are surprisingly similar, albeit quantitatively larger interaction strengths 
are required in three dimensions for getting a Mott-Hubbard insulating state. The exchange pa- 
rameters of an effective Kugel-Khomskii-like spin-orbital-model are also derived and reveal strong 
antiferromagnetic tendencies. 
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I. INTRODUCTION 



With the tremendous experimental progress^ in grow- 
ing oxide heterostructures, new routes to engineer the 
dimensionality of materials have become available. A 
physical phenomenon particularly sensible to dimension- 
ality is high-Tc superconductivity^ which seems to find 
its ideal playground in quasi two-dimensional systems. 
Maybe not accidentally, in two dimensions also strong 
antiferromagnetic spin fluctuations prosper, without un- 
dergoing a phase transition to long-range order as in three 
dimensions. 

Besides two dimensionality, the electronic structure of 
the most prominent class of unconventional superconduc- 
tors, the cuprate family, is governed by a single Cu (1^2 _y2 
orbital doped slightly away from half-filling. An impor- 
tant step towards a better understanding of cuprate su- 
perconductivity was taken by Pavarini et al.^ who found 
an empirical trend in the material dependence of the crit- 
ical temperature, Tc, at optimal doping. In the proposed 
scenario the material dependent parameter is the energy 
difference between the planar Cu d^2_y2 orbital at the 
Fermi energy and an axial Cu 4s-like orbital at higher 
energy. 

In a recent Letter'^, Chaloupka and Khaliullin pro- 
posed the artificial engineering of nickelate heterostruc- 
tures with the aim of getting cuprate-like electronic struc- 
ture and interplay between axial and planar orbitals. A 
possible realization are LaNi03=LaO-Ni02 layers alter- 
nated with insulating layers such as LaA103=LaO-A102 
with the same sequence of sub-layer charges. To avoid 
macroscopic polarization, the substrate on which the het- 
erostructure is grown should preferably have the same se- 
quence of sub-layer charges (SrTiOs may for instance not 
be suitable) . We have shown that in these d^ heterostruc- 
tures, the Ni (i3j,2_j.2 orbital can play the role of the axial 
orbital^; it is pushed considerably above the planar Ni 
dT.2_,,2 orbital as an effect of either the Coulomb inter- 



action, an axially compressive strain, or control of the 
apical cation (substitution of Al by, e.g., Sc)^ so that a 
single Fermi surface sheet with similar geometry and vol- 
ume as in cuprates emerges, before the system becomes 
eventually insulating. 

The main focus of this complementary paper is a de- 
tailed analysis of the Mott-Hubbard transition taking 
place in the axial /planar two-band system and of the 
changes observed when going from the two dimensional 
heterostructure to three dimensional bulk nickelates. 
The latter can be viewed as a dimensional crossover, in 
which the Cu d^^2_j.2 degrees of freedom are more equiva- 
lent to the Cu dx2-y2 ones due to the lack of dimensional 
constraint. 

For the purpose of this study, we simplify the local- 
density-approximation (LDA) bandstructure of Ni eg or- 
bitals by including merely the hopping between the first- 
and second-nearest neighbor Ni d(eg) Wannier orbitals, 
in contrast to Refj^ where the full LDA bandstructure 
was considered. This gives a quarter-filled two-band 
model for the eg degrees of freedom, i.e., one electron 
in the two Ni eg orbitals. We recall here that in transi- 
tion metal oxides the transition metal ion is usually sur- 
rounded by an oxygen octahedron and, depending on the 
specific case, either t2g (as, e.g., in vanadiumsesquioxide) 
states or eg states (as in the nickelates and cuprates) 
play the most important role in determining the low 
energy excitations of the system. Hence, our study of 
an Cg-band model at quarter filling and the dimensional 
crossover from quasi-two- to three-dimensional systems 
is of interest also from a more general perspective. More- 
over, since there are plenty of parameters governing the 
physics such as the bandwidth of the two bands, their 
crystal field splitting, the interband hybridization, the 
Coulomb interaction, and last but not least the dimen- 
sionality, it is of great importance to derive such models 
from realistic bandstructures. 

This way we also make some contact to previous multi- 
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band model studies^"— where features such as the orbital 
selective Mott transitionsi^"— were analyzed. Many of 
these multi-band studies made approximations, such as 
neglecting the hybridization between the orbitals. This 
poses the question in how far the employed models are 
applicable to real world materials. Moreover, most of the 
aforementioned studies concentrated on half filling. 

In Section |TT1 we will briefiy discuss the two-band 
model and how we derive its parameters. In Section Hi Bl 
we discuss the evolution of the spectral function, self en- 
ergy and double occupation with increasing Coulomb in- 
teraction strength. In Section [llll we turn to the dimen- 
sional crossover, comparing the quasi two-dimensional 
system for a heterostructure and a more homogeneous 
3D model. Finally, we discuss the issue of magnetic order 
at low temperatures and derive an effective spin-orbital 
superexchange model in Section IIVI before we arrive at 
the conclusion in Section IVl 



II. LDA+DMFT 
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FIG. 1. (Color online) Orbitally projected DOS (spectral den- 
sity) for each of the two eg bands (thick lines). Neglecting 
the inter-orbital hybridization leads to the results shown as 
thin hnes. The hybridization is seen to strongly enhance the 
width of the Zz^ — spectral density, leading to a remark- 
able itinerancy, despite the quasi two-dimensional nature of 
the LaNiOa /LaAlOa heterostructure. 



In the basis {\x'^ — y'^), \3z'^ — r^)} of the planar and 
axial Cg Bloch-summed orbitals, the one-electron part of 
the Hamiltonian is: 



A. Effective two-band model 

Density functional calculations were performed 
for the simplest, 1/1 unstrained super lattice 
LaNiOa/LaAlOa^LaO-NiOz-LaO-AlOa of which a 
detailed discussion can be found in Rcf. 5. Here we 
recall the most important facts relevant for the dis- 
cussion of our model study. In the LaNiOs/LaAlOs 
heterostructure nickel is in a threevalent state Ni'^"'" . The 
coordination polyhedron is an octahedron with a small 
tetragonal distortion due to the LaAlOa layers along the 
c-axis. The cubic crystal field eigenstates of the Nickel 
d-electrons are the lower lying t2g triplet and the higher 
lying Cg doublet. The cubic splitting is large enough 
so that the nominal filling in Ni'^^ is (completely 
filled) and (quarter filling). As a consequence we 
find two eg bands at quarter filling around the Fermi 
energy to be the relevant degrees of freedom for the 
low energy excitations of the system^. The insulating 
LaAlOs layers give rise to the tetragonal distortion 
which splits the axial |3z^ — r^) orbital ~ ISOmeV above 



the 



y ) orbital and, more importantly, strongly 



reduce the width of the axial \3z'^ — r^) band by blocking 
the perpendicular hopping^. 

In order to obtain the model for our DMFT study we 
construct a basis of localized Ni Cg Wannier functions 
using the Nth-order muffin-tin orbital method (NMTO) 
like in Ref. but now truncate the hoppings after those 
between second-nearest Ni neighbors. For our purposes 
of model calculations we prefer this simpler picture of 
the systems kinetics in real space and pay the price of 
a slightly less accurate representation of the dispersion 
of the LDA conduction bands. Let us remark that this 
simplification could affect the Fermi surface transition 
discussed in Ref. 5 but not the overall properties of the 
MIT and the 2D to 3D crossover. 
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Here, energies are in eV. The first and second terms ac- 
count for respectively the intra- and inter-orbital nearest- 
neighbor hoppings. We see that the hopping between the 

— y^ orbitals is almost three times as large as the one 



between the 32 



orbitals and that the inter-orbital 



hopping is nearly the geometrical average of the intra- 
orbital hoppings. This, as well as the factor 3, is what 
one expects when the Ni-Ni hopping is exclusively via 
O, i.e., when iNi-Ni = ^Ni-o • to-Ni/(eNi - ^o) and tm-o 
follows the canonical rules^^, which in the present case 
yield: t^-i^y-i — ^/St^z^-r^ in a;-, and y-direction. The 
third term in equation ^ is the second-nearest neigh- 
bor hopping, and the fourth is the on-site term, i.e. it 
gives the tetragonal crystal-field splitting. The two last- 
mentioned terms do not mix the two orbitals. In real 
space, the Hamiltonian of our effective model reads: 



mod. 



int. ; 



(2) 



where l,m d {p, a} label the planar and axial orbitals. 
The hopping amplitudes tnjm are the matrices in Eq. ([1]) 
with l,m as the orbital index and i,j as site index of 
nearest or next nearest neighbors on the square lattice. 

The effects of the inter-orbital hybridization are clearly 
visible in Fig. [T] showing the noninteracting density 
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of states resulting from the noninteracting part (jlj of 
Hamiltonian A comparison to the case without hy- 
bridization (tpa = 0) plotted as a thin line shows a big 
difference especially for the 3z^ — r^-projection. Its width 
is strongly enhanced if the hopping to the more mobile 
x^— y^-orbitals is properly taken into account. This effect 
is important since otherwise the reduced hopping in the 
z-direction would more severely reduce the itinerancy of 
electrons in the 3z^ — r^ orbitali^. That is, in the absence 
of the x'^ — y^-to-Sz^ — hybridization the 3z^ — or- 
bitals would localize too easily (even at a small Coulomb 
interaction strength). This scenario was originally pro- 
posed in Ref.'^ prior to realistic material calculations. 
The Coulomb interaction terms read: 

H,,,t.^Uj2nlni+ J2 {V-5.,.'J)n1inl'^ (3) 

il(T il^nicrcr' 

where nfj — c\i^Ciia, So-, a' denotes the Kronecker sym- 
bol, U and V represent the intra-orbital and inter-orbital 
repulsion, respectively, and J is the Hund exchange term. 

For the solution of the impurity problem in the DMFT 
self consistency loop we use quantum Monte Carlo 
(QMC) simulations of the Hirsch and Eye typei^ where 
neither the pair hopping term (corresponding to highly 
excited states) nor the spin-flip term (which poses a se- 
vere sign problem in the numerical simulations^,) have 
been included. Presented results were obtained at /3 — 
10eV~^ except for the double occupancies shown in the 
right panel of Fig. [5j Further, convergence with respect 
to the value of the Trotter discretization At was vali- 
dated. Afterwards, the maximum entropy method was 
employed to perform analytic continuation of the Green 
function to the real axis. Concerning the choice of the 
interaction parameter we perform a systematic analysis 
where U is varied in a regime of physically meaningful 
values in order to investigate the Mott-transition and 
its specific features. The Hund's exchange term, on the 
other hand, has been fixed to J = 0.7eV, a reasonable 
value close to that of atomic Ni. Moreover, assuming the 
deviation from cubic symmetry to be negligible with re- 
spect to the Coulomb interaction parameters, we employ 
V = U-2J. 

Before presenting our numerical results a comment 
about the applicability of DMFT to the finite dimen- 
sional system considered is due. While it has been shown 
that corrections arising from spatial correlations beyond 
DMFT— are generally small in 3D (except in the close 
proximity of a second order phase transition), stronger 
effects have to be expected for 2D systems. The strength 
of such correction with respect to DMFT is, however, de- 
creasing with increasing temperatures^SiS^ and therefore 
can be neglected as a first approximation in the temper- 
ature regime we are considering here. This corresponds 
to the assumption that we are well above any ordering 
temperature at /3 = 10 — 25eV~"'^. Naturally a differ- 
ent treatment to include effects of spatial correlations is 
necessary in order to study the low temperature regime 
where spin fiuctuations and, possibly, superconductivity 



emerge. For this reason we include in Sec. IIVI a two-site 
calculation to estimate the size of non-local antiferromag- 
netic correlations for the system considered. 

B. Mott transition in two dimensions 

Spectral Function and Self Energy 

Fig. [2] shows the evolution of the fc-integrated spectral 
functions A{uj) for both orbitals with increasing U. Let 
us first focus on the x'^ — y^-projected spectral density 
(black curve). Already at f7 = 4.4 eV {V = 3.0 eV; 
first panel), the spectrum is strongly renormalized with 
respect to the noninteracting case: A narrow quasiparti- 
cle peak at ep and Hubbard bands which are separated 
by an energy of the order of U are well visible. As the 
value of the interaction parameter is increased to U — 5.4 
eV and U = 6.4 eV (second and third panels), such 
spectral features become even more pronounced as one 
observes a stronger renormalization of the quasiparticle 
peak and more pronounced Hubbard bands. However, 
a,t U — 7.4 eV (fourth panel) a qualitative change has 
occurred: the spectral weight at the Fermi level has van- 
ished and a Mott-Hubbard gap has been formed. Hence, 
a Mott metal-insulator transition has taken place be- 
tween U = 6.4 and 7.4 eV. 

Fig. [2] also indicates that the metal-insulator transi- 
tion of the 3z^ — r^-band (red curve) occurs for the same 
critical interaction oi U ^ 6.4 eV as in the x'^ — band. 
It is, however, qualitatively different. For V increasing 
from 4.4 eV, one observes that the 3z^ — quasiparticle 
peak is shifted above that of the x^ — j/^ peak, which is 
at the Fermi level. This is the Coulomb enhancement of 
the tetragonal crystal field which causes the Fermi sur- 
face to have only one, mostly x^ — j/^-like sheet as in the 
cuprates^'. For increasing U, the 3z^ — spectral weight 
shifts further and further above the Fermi level where for 
U — 1 A eV it has completely disappeared. 

The remaining spectral weight below the Fermi energy, 
corresponding to the lower Hubbard band, is entirely due 
to the hybridization with the lower x^ — y^-like Hubbard 
band. 

The different nature of the metal-insulator transition 
in the x^ — y^- and the 3z^ — -bands can be identi- 
fied even more clearly from an analysis of the self en- 
ergy. In Fig. [3j we show the evolution of the imaginary 
part of the self energy ImE3,2_y2 (a;„) with increasing in- 
teraction. Again, for the x^ — band, the expected 
behavior of a standard Mott transition is found: The 
appearance of the insulating gap results from the diver- 
gence of Iml]2,2_y2 (w) in the limit of zero frequency. In 
this limit the quasiparticle renormalization factor Z ~ 
(1 — ^^^^^7^ |w=o)~^ goes to zero as a further hallmark 
of the Mott character of the transition in the x^ — y^- 
band (see black circles in the left panel of Fig|31). At 
the same time, however, the insulating transition of the 
3z^ — r^-band is not associated to any qualitative change 
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FIG. 2. (Color online) Evolution of the LDA+DMFT spec- 
trum with increasing Coulomb interaction U (with J — 0.7eV, 
V = U - 2J) mt 13 = lOeV"^ The black curves show the 



X ~y and the red/gray ones the Zz 
density. 



- r -projected spectral 
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FIG. 4. (Color online) Left hand side: Effective chemical 
potential ^ — ReE(O) relative to the noninteracting chemical 
potential for the —y^ (black) and —r^ (red/gray) orbital 
at /3 = 10eV~^. At the Mott transition, the low frequency 
effective chemical potential of the 3z^ — orbital is strongly 
reduced so that it is shifted above the Fermi energy, i.e., the 
3z^ —r^ orbital becomes "band" -insulating. Right hand side: 
Quasiparticle weight Z for the — (black) and 3z^ — 
(red) orbital (also at /3 = 10eV~^). For the x'^ —y^ orbital, we 
find Z — > 0, i.e., the (central) quasiparticle peak disappears 
and the x^ — orbital becomes Mott-insulating. 
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FIG. 3. (Color online) Evolution of the LDA-^DMET self 
energy (vs. negative Matsubara frequencies) with increasing 
[/ at /3 = 10eV~^. The black curves show the x^ — y'^ and 
the red/gray ones the — self energy. Via extrapolation 
of the imaginary and real part of the self energy to a;n — >■ 
we can obtain the quasiparticle renormalization Z and the 
"crystal-field" enhancement, respectively, shown in Fig. [l] 



in ImS322_r2 (w). For all interactions a "metallic"-like 
bending can be observed in ImSi-^^i (lo) at low frequen- 
cies. As a consequence of this, the renormalization fac- 
tor Z stays always finite (red crosses in the right panel of 
FigH]) . The disappearance of spectral weight at the Fermi 
energy in this case is determined by a large relative shift 
of the orbitals induced by the Coulomb interaction. This 
shift A^, is a consequence of the strong variation of 
the real part of the self energies at zero frequency shown 
in the left panel of Fig. H) More specifically, while the 
effective chemical potential of the — y^-band (defined 
as — ReE^2_j,2 (0)) remains close to the noninteracting 
value up to [/ = 6.4 eV, its value for the Sz^ r^-band 
is clearly decreasing with increasing U . As a result, the 
metal insulator transition of the — r^-band occurs 
for a value oiU ^ 6.4 eV, when the low energy spectral 
weight of this band is shifted above the Fermi energy. 

A key feature in the physics we have observed here is 
the interplay between the two orbitals stemming from the 
fairly large hybridization in equation proportional to 
cos kx — cos ky . An evident sign of such interplay is the 
concomitance of the metal insulator transition in both 
orbitals. For understanding the underlying physics it is 
instructive to consider the different possible scenarios for 
the MIT in a two band system with hybridization. 

Increasing the Coulomb interaction has basically two 
effects (see also Refs. HI, HI and in particular Ref. IT^ : 
Double occupancy of a site becomes more expensive and 
on-site energies of the two orbitals are shifted with re- 
spect to each other. This gives rise to three possible 
transition scenarios, (i) If the shift between the bands is 
increasing slowly with interaction, double occupancy be- 
comes practically forbidden for a value of V where both 
bands are still partially filled. As a result, a simultane- 
ous Mott transition occurs with Z — for both bands. If, 
instead, the shift is increasing more rapidly with interac- 
tion, at a certain point one of the bands will empty, leav- 
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ing the other one half filled. The nature of the transition 
then depends on the value of U for which this happens. 
There are two possibilities: (ii) If U is strong enough 
to forbid double occupancy in the half-filled band, this 
band will become insulating abruptly. Such a situation 
is reflected in simultaneous, but qualitatively different, 
metal-insulator transitions for both bands. This is what 
we observe in Figs. [5] and [3] (iii) The other possibility 
is that U is not yet strong enough to prevent the double 
occupancies in the half-filled band. The Mott transition 
in this band - differently from our case - would then 
not take place simultaneously with the depletion of the 
other band. This case is different from the one considered 
in Figs. [2] and [3] and the only one in which the metal- 
insulator-transition can be fully described in terms of an 
effective single-band model. 

The initial lifting of the orbital degeneracy (the LDA 
crystal-field splitting) is expected to push the system in 
the direction of situation (ii) , which we observed for our 
model with an initial crystal field splitting of £3^2 — 
ex:2_y2 = 150 meV. 

The inter-orbital hybridization also gives rise to two ef- 
fects. The first one, which can be understood intuitively, 
is an effective broadening of the more narrow band (see 
Fig. [1] for the noninteracting case). This leads to a more 
metallic behavior of the — r^-electrons. The second 
effect is more intrinsic in the sense that the hybridiza- 
tion makes the bands more similar, obviously pushing 
the system towards situation (i). It is important to no- 
tice that both effects work against the effective crystal- 
field splitting. As mentioned above, the fingerprints of 
the hybridization can clearly be seen also in Fig. [21 Al- 
though the narrow 82^ — -like-band is shifted above 
the Fermi energy, at the metal-insulator transition some 
residual 3z^ — weight remains in the region of the lower 

— y^-like-Hubbard band. 



Double occupancy 

As in the case of the half-filled Hubbard model, com- 
plementary information about the metal-insulator tran- 
sition can be extracted from analysis of the double occu- 
pancy d, which is the derivative of the free energy with re- 
spect to the interaction parameter U and is a good indica- 
tor for the Mott metal-insulator transitior^. However, 
such an analysis for a two-band model is more compli- 
cated since it involves not only the orbital-diagonal dou- 
ble occupations (dmm = (?^mt'^m^)) but also the orbital- 
offdiagonal parts with parallel (djj^ = (np^naf)) and an- 
tiparallel spin orientations (d^^ = [np^Ua]) )■ 

The evolution of these four quantities is shown in Fig. [5] 
for two different temperatures. In both cases we observe 
the same trend for all double occupations which are de- 
creasing with increasing interaction. Starting the anal- 
ysis with the — j/^-orbital, which undergoes a Mott- 
Hubbard transition when getting half filled, we observe 
that its double occupation decreases around U — 6.4 eV, 
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FIG. 5. (Color online) Double occupancy vs. Coulomb inter- 
action U for two values of P (left: lOeV"^; right: 25eV"^). 
Shown are the double occupation of the — orbital (dpp), 
of the — orbital (daa) as well as the double occupa- 
tion of both orbitals with parallel (dJJ) and antiparallel spin 
(dp„). Inset: double occupancies normalized by the respective 
density {Smn see text). 



i.e., near the critical value of the interaction that we 
can estimate form the spectral function in Fig. [21 At 
kT= 1/lOeV (Fig. [SI left panel), however, this behavior 
is thermally smeared out. The Mott-Hubbard transition 
occurs in form of a very smooth crossover. A better es- 
timate can be made by going down in temperature to 
kT= l/25eV (Fig. [51 right panel), where we can see a 
very steep drop of dpp marking the transition point at 



for the 



y 



U sa 7.0 eV where we also found Z 
band. 

At a first glance the behavior of the other three double 
occupancies daa, d^i, and appear to be more diffi- 
cult to interpret: One could have expected that djl and 
d^^ should be the largest ones because they are energet- 
ically less expensive {V^V — J) in comparison with the 
dmm {U). What we observe, however, is a crossing of dJJ 
with dpp 'At U Ri 5.4 eV and U « 4.2 eV for ;9 = 10 and 
f3 ~ 25eV respectively. Moreover, dj,^ is lower than dpp 
for all interaction values. This situation can be properly 
understood considering also the depletion of the 3z^ — ri- 
band: The low values of daa, dpi, and djjj reflect the low 
3_j,2 _ 7-2_giec^]-on density. In order to disentangle these 
effects we have plotted the double occupancies normal- 
ized by the respective density in the inset of the left panel 

of Fig. MiSmni = dmm/(?^m)^ = d^H {{Up^) {Uai)) , 

and S^a — dpX/{{^pf){''^af))- In this way the hierarchy 
of the energetic consideration is restored and all renor- 
malized double occupancies behave in a similar way. On 
the other hand, the crossing observed in the unrenormal- 
ized data can be also interpreted as a clear signal that 
we are far from the situation (i) described above, that is, 
a simultaneous Mott-Hubbard transition in both bands. 



III. FROM 2D TO 3D 

Having discussed in detail the results for the two- 
dimensional model for the LaAlOs/LaNiOa superlattice, 
we now present results for a 3D model for bulk LaNiOa, 
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FIG. 6. (Color online) Same noninteracting projected DOS as 
Fig. [S] but for the 3D case. The — band is considerably 
broader than in the layered 2D case due to the additional 
hopping in the z-direction. 



and compare with the 2D case. Specifically, we no longer 
restrict the hopping of the axial Zz^ — r^-orbital to the 
xy-plane, but allow for hopping along the z~axis by 
adding the term 



-2 -tz ■ cos{kz 




1 



to the kinetic part of our Hamiltonian ([Ij. We choose 
tz ~ 0.6eV, i.e., 4/3 times the 0.45 eV in equation (P), as 
it should be in cubic symmetry according to the Slater- 
Koster tables^! and Ref. [H. Our 3D model is not com- 
pletely cubic because from equation (jlj wc kept the 
tetragonal crystal-field splitting (last term in Eq. ([1])) 
as well as the second-nearest neighbor hoppings in xy- 
plane, without including them in the xz- and yz-planes. 
Moreover, this 3D model is not to be seen as a true 
LDA+DMFT analysis of pseudo-cubic LaNiOa, albeit 
the link to this compound is a motivation for the analy- 
sis of the modified model. 

The noninteracting density of states for this 3D model 
is shown in Fig. [S] The main difference from the 2D 
model in Fig. [1] is obviously, a much broader 3z^ — ri- 
band caused by the additional hopping in the 3rd direc- 
tion. As a result of this broadening the two different 
bands are much more similar than in the 2D case. In the 
truly cubic case the x^ — y^ and 3z^ — spectral densities 
would have been identical. 

In Fig. [7] we show the spectral functions and self ener- 
gies of the interacting 3D system. On the left hand (right 
hand) side the data for a metallic (insulating) solution are 
plotted. They refer XoU — lA and U = 9.9 eV, respec- 
tively. This indicates that the metal-insulator transition 
occurs for considerably larger values of the interaction 
than in the 2D system, as might be expected. Another 
difference to the 2D case is the enhanced spectral weight 



of the iz'^ 



-band below the Fermi energy. However, 



the metal-insulator transitions for the two bands are 
qualitatively similar to those of the 2D case as is clearly 
demonstrated by the behavior of the imaginary part of 
the selfenergy shown in the bottom panels of Fig. [T] 
The evolution of the self energies qualitatively resembles 
the one shown in Fig. [3] displaying the same "metallic" 
down-bending. The only (quantitative) difference to be 
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FIG. 7. (Color online) Evolution of the spectrum and the 
LDA-fDMFT self energy (vs. Matsubara frequencies) with 
increasing Coulomb interaction V = \J — 2J for the 3D case 
at /3 = lOeV"^ The black curves show the ■3? — spectrum 
and the red/gray ones the ?>z^ — spectrum. 



noticed w.r.t. the 2D case is a larger imaginary part 



of the 3^2 



-selfenergy at low-frequencies. This can 



be understood as a natural consequence of the closer re- 
semblance of the bands in the noninteracting density of 
states. In terms of our precedent discussion we thus ob- 
serve again the situation (ii) of section fllBi i.e., a simul- 
taneous albeit quantitatively different Mott transition. 
However, we are now slightly closer to the case (i) of a 
simultaneous double Mott transition for both bands than 
for the layered model. In other words, the Mott transi- 
tion is qualitatively the same for the 3D as the 2D case, 
but quantitatively the Mott transition is shifted to much 
higher values of the interaction strength. Such large val- 
ues can hardly be realized in experiment. 

It must however be emphasized that for a truly cu- 
bic 3D model, the two Cg orbitals remain identical (for 
all strengths of the interaction) so that the Mott transi- 
tion will occur simultaneously for both bands and for a 
value of U which is about y/2 larger than for the single- 
band easels. That a small symmetry breaking at the 
LDA level, as exhibited by our 3D model, can dras- 
tically change the nature of the Mott transition in the 
paramagnetic state has previously been found in DMFT 
calculations for 3D systems2ii2^ and t\ systems'^", as 



well as for the e\ system LaMnOg^. 



IV. PREDOMINATING SPIN FLUCTUATIONS 
AND EFFECTIVE SUPEREXCHANGE MODEL 

Predominating spin fluctuations 

The model calculations discussed in the previous sec- 
tions were performed for the paramagnetic phase. How- 
ever, one expects the system to display magnetic and/or 
orbital order at sufficiently low temperatures. The treat- 
ment of such complex ordered phases is particularly dif- 
ficult in LDA-I-DMFT. An insight into this physics can 
be obtained however by diagonalizing a "two-site ver- 
sion" of Hamiltonian (with open boundary conditions, 
i.e., a diatomic molecule directed along the a;-axis). The 
analysis of the corresponding energy levels for this two- 
electron system allows to infer the relevant fluctuations 
which dominate the low temperature physics. On the 
left-hand side of Fig. [8] we show the calculated energy 
level diagrams as a function of U for the same set of 
parameters that we have used for the LDA-I-DMFT 2D 
crystal calculation discussed in previous sections. 

For the values of U considered, the mixing of the low- 
est energy levels with the energy-costly double-occupied 
states is very small, i.e., these states are already above 
the energy window which we plot in Fig. |51 More inter- 
estingly, our results show that the ground state of the two 
site system is always a spin singlet with the two electrons 
residing predominantly in the same orbital (a:^ — y^) on 
each site. The splitting of the singlet ground state and 
the excited triplet is large, and, in the large U limit, it 
can be regarded as an exchange integral of Jcx. ~ 200 
meV. This indicates an instability of the system at low 
temperatures towards a spin antiferro- and orbital ferro- 
ordered state. We have also studied the effect of varying 
the initial level (i.e., crystal field) splitting for a reason- 
able value oiU — 5eV (right-hand side of Fig. IS]). The 
main effect of an enhancement of the crystal field splitting 
is to strengthen the tendency towards ferrotype orbital 
ordering (of the — and the — states) even- 
tually leading to purely one-band physics, with only the 



^ orbital involved. 



y 



Effective superexchange model 

A first analysis of the magnetic and orbital ordered 
phases in the La2NiA106 heterostructure was carried out 
in Ref. 4 considering the following superexchange Hamil- 
tonian, defined for a bond ij \\ 7 in the Ni02 plane: 



R. 



CT 
■a,±,± 



)( 



,1 



(7) 



)P..^J (4) 



with an implied sum over ct = 0, 1 and all combinations 
of ±, ±. Pa,ij represents the projector to a singlet and 
a triplet state of two Ni'^^ S = ^ spins while ^ ± fj''-' 
selects the planar orbital in the plane perpendicular to 
the 7 axis and the directional orbital along this axis. 



> 

" -0.2 




2.5 3.5 4.5 5.5 6.5 7.5 0.0 0.1 0.2 0.3 0.4 0.5 

U [eV] AcF [eV] 

FIG. 8. (Color online) Left hand side: Energy level diagram 
of the two site model as a function of the Coulomb repulsion 
U for AcF = 0.15eV (dashed line in the plot on the right- 
hand side) . Right hand side (reproduced from Ref. @) : Same 
diagram as a function of Acf for U = 4.5eV (dashed line 
in the plot on the left-hand side). The color of the points 
indicates the character of orbital ordering (of: orbital ferro; 
Oaf: orbital antiferro) 



Neglecting the charge transfer part of Hamiltonian 
allows for a comparison with our two-site calculation, al- 
beit only for Acf = OeV since otherwise Hamiltonian 
(HI has another orbital symmetry than the full Hamil- 
tonian. In this way we can perform the first numerical 
estimate of the Ra.±.± coefficients of the superexchange 
Hamiltonian Q as a function of the parameter of the mi- 
croscopic Hamiltonian. More specifically, the Rcr,±,± co- 
efficients of Fig. ini are calculated form the level-splitting 
seen in Fig.[8]which, in turn, depends on the values of the 
parameters U, J, and t of Hamiltonian ©. The result 
of this comparison are shown in the left panel of Fig. ^ 
where for simplicity we assumed the two orbitals to be 
degenerate. 

It is important to note here, that the orbital projectors 
in Eq. do not refer to the — y'^ and 3z^ — basis 
we have used so far since the 7 axis is lying in the Ni02 
plane. This is also refiected in the mixed color of the 
R°__ coefficient which corresponds to the singlet ground 
state. 

Finally, in the right panel we show the effect of turn- 
ing on the crystal field splitting Acf, disregarding some 
(here still small) terms which admix the orbitals in 
Hamiltonian (U). The changes in the color clearly indi- 
cates a rapid demixing of the eigenstates, since the x^—y^ 
and 3z^ — r^ orbitals are eigenfunctions of the crystal field 
operator. This limits the applicability of Hamiltonian 
to the cases of small crystal field splittings. 

As for the magnetic fiuctuations, the Kugel-Khomskii- 
type Hamiltonian clearly favors an antiferromag- 
netic spin-alignment along each pair. Depending on the 
crystal-field splitting, the preferred orbital orientation is 
along the bond or -for larger crystal field- a ferro-orbital 
occupation of the — y^ orbital. 
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FIG. 9. (Color online) Left hand side: Evaluation of the R- 
coefficients of Hamiltonian Q as a function of the interaction 
parameter U calculated with the two site model in the — 
and 3z^ — basis for Acf = OeV . Right hand side: Evolution 
with increasing crystal field splitting starting from the largest 
U value of the left panel {U — 7.5eV). The color of the points 
indicates the character of orbital ordering with the same scale 
as in Fig. \8\ 



V. CONCLUSION 

We have analyzed in detail the Mott transi- 
tion in the high-temperature, paramagnetic state of 
LaAlOa/LaNiOs heterostructures, including the dimen- 
sional crossover from the two-dimensional layered het- 
erostructure to the three dimensional bulk system. The 
LDA bandstructure was NMTO-downfoldcd to the Ni Cg 
orbitals and, in order to concentrate on the essentials of 
the Mott physics, we included merely the hoppings be- 
tween first and second-nearest Ni neighbors. This realis- 
tic model also allowed us to contribute to the intensive 
discussion on the orbital-selective Mott transitior^^rii, a 
discussion which with a few exceptions^ has hitherto 



neglected the effect of interband-hybridization. This hy- 
bridization turned out to be most important for the itin- 
erancy of the 3z^ — r^ band and for a simultaneous, albeit 
different Mott transition in both orbitals. 

For realistic U values we find that LaNiOa/LaAlOs is 
a metal with an (a;^ — ?/^)-like quasiparticle peak at the 
Fermi level and a {3z'^ — i/^)-like peak slightly above. 
At larger values of the Coulomb interaction, the 3z^ — 
spectral density is shifted above the Fermi as in a metal to 
band-insulator transition for this band, while the x^—y^- 
like band undergoes a genuine Mott transition with diver- 
gent effective mass (XjZ). Due to the hybridization be- 
tween both orbitals, there is however still spectral weight 
in the Sz^ — channel at the position of the predomi- 
nately x^ — y^ lower Hubbard band. 

If LaAlOa/LaNiOs is grown on a substrate, stress can 
further enhance the initial (LDA) orbital splitting so that 
one could obtain a situation which is more single-band- 
like. This is expected to reduce the critical Coulomb 
interaction strength for the Mott transition. 

Studying the two-site version of the model, which also 
includes second order processes in the nearest-neighbor 
hopping, we made contact with an effective Kugel- 
Khomskii-Hamiltonian, indicating predominant antifer- 
romagnetic spin fluctuations, an essential ingredient for 
a prospective superconductor. Antiferromagnetic long- 
range order is hence also possible and might enhance the 
tendencies towards an insulating ground state. 
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